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Abstract. 

Flow coefficients v n for n = 2, 3, 4, characterizing the anisotropic collective flow in 
Au+Au collisions at ^/snn — 200 GeV, are presented. They indicate the expected 
growth of viscous damping for sound propagation in the quark gluon plasma (QGP) 
produced in these collisions. Hydrodynamical model comparisons which include the 
effects of initial state geometry fluctuations, highlight the role of higher harmonics 
{v n ,n>2) as a constraint for disentangling the effects of viscosity and initial conditions, 
and suggest a small specific viscosity for the QGP. This viscosity is compatible 
with that obtained via a newly proposed technique [6] which employs the relative 
magnitudes of v n to estimate the viscosity, and the "viscous horizon" or length-scale 
which characterizes the highest harmonic that survives viscous damping. 

1. Introduction 

Measurements of anisotropic flow in heavy-ion collisions at the Relativistic Heavy Ion 
collider (RHIC), continue to play a central role in ongoing efforts to characterize the 
transport properties of the quark gluon plasma (QGP) produced in these collisions. 
Recently, considerable attention has been given to extractions of the specific viscosity 
r]/s [the ratio of viscosity (77) to entropy density (s)] via elliptic flow (V2) measurements. 
These extractions indicate a range of 1-2 times the conjectured lower bound [2] for rj/s, 
(i.e. ~ 1 — 2). This large uncertainty (100%) is known to be dominated by the 
uncertainty in model estimates of the initial eccentricity [3j H]. Thus, a more precise 
extraction of r]/s requires the development of new experimental constraints. 

Because of the acoustic nature of flow (i.e. it is driven by pressure gradients), 
a transparent way to evaluate the strength of dissipative effects, is to consider the 
attenuation of sub- horizon sound modes in the plasma. In the presence of viscosity (77), 
sound intensity is exponentially damped as : 5T^ v (t) = exp (— |||^) (57^(0) [5], where 
the spectrum of initial (t = 0) perturbations of the energy- momentum tensor T^, can be 
associated with the harmonics of the shape deformations and density fluctuations of the 
collision zone; k is the wave number for these harmonics, and t and T are the expansion 
time and the temperature of the plasma respectively. Since viscous damping scales as 
k 2 the viscous corrections for the eccentricity driven harmonics v n (with wavelengths 
2iiR/n for n > 1, i.e. k ~ n/R), are expected to scale as n 2 K; R is the transverse 
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size of the collision zone and K is the Knudsen number [BJ. The latter is often used 
to parametrize viscous corrections. The length scale r v or "viscous horizon" separates 
the sound wavelengths which are effectively damped out, from those which are not, and 
k v = 2ir/r v is linked to the order (n v ) of the highest harmonic which survives viscous 
damping. Thus, the relative magnitudes of the higher-order harmonics (i> n ,n>3) are 
expected to provide additional constraints on both the magnitude of rj/s and the "best" 
model for eccentricity determinations [5J [7J [10] . Here, we report new v n measurements 
PQ and investigate their utility as constraints for precision extraction of rj/s. 

2. Data analysis and results 

The results presented are derived from ~ 3.0 x 10 9 Au+Au events obtained during 
the 2007 RHIC running period. The Fourier coefficients v n were obtained via (i) 
pair-wise distributions in the azimuthal angle difference (A0 = <ft a — <pb) between 
particles with pseudorapidity separation At/ > 1 and transverse momenta and 
p h T (respectively); dN J^ s oc (l + J2n=i ^ v n(PT) v n(PT) cos(nA0)) , and (ii) azimuthal 
distributions oc (1 + J2n=i 2v n cos(n</> — n^ n )) , of charged hadrons detected in the 
PHENIX central arms at an azimuthal angle 0, relative to event planes ty n [H] obtained 
with three separate detectors: Beam-Beam Counters (BBC), React ion- Plane Detectors 
(RXN), and Muon Piston Calorimeters (MPC). Each detector has a North (South) 
component which allows correlation studies between sub-event planes $ n determined at 
forward (backward) rapidity. The absolute pseudorapidity coverage for these detectors 
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Figure 1. Correlation strengths of the event planes for various detector combinations 
as a function of collision centrality. The detectors in which the event plane is measured 
are: RXN North (A), BBC South (B), MPC North (C) and MPC South (D). 

Figure [TJ shows the centrality dependence of the correlation strength 
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^cos(j($^ — f° r sub-event combinations (A, B) involving different event-plane 

detectors with Arj' ~ 5 and Arf ~ 7. Note that the order j is chosen to account for 
the n-multiplet of directions (2?r/n) of and the magnitudes for ^sin(j($^ — 
are consistent with zero for all centrality, j, and $ combinations. Positive sub-event 
correlations are indicated in panels (a) and (b) for ^2,3,4, with sizable magnitudes for 
^2,3 and much smaller values for ^4. The negative correlation indicated in panel (a) is 
due to the well known antisymmetric pseudorapidity dependence of sidewards flow (t>i), 
as well as momentum conservation. The expected correlation between ^ and ^2, and 
that between ^2 and \l/4 are confirmed in panel (c); they show that ^2 and ^4 are 
correlated with the reaction plane. An initial state fluctuation origin of ^3 [and hence 
1*3] is well supported by the absence of a correlation between \l>2 and ^3 in panel (d). 
The absence of a correlation between fy 2 and ^3 reflects the rather large fluctuations of 
\l/3 about ^2 and gives a null value for v 3 measured relative to |H]- 




Figure 2. v n vs. pr for several centrality bins as indicated. 



Figure |2] shows a clear growth of V2 from central to mid-peripheral collisions, and 
a near constancy of V3 and V4 which are strong indicators of the role of initial state 
fluctuations in establishing the higher harmonics. It is noteworthy that the ratios 
f 3 /(f 2 ) 3//2 and v 4: /(v2) 2 are essentially independent of pr (for p T < 2.5 — 3.0 GeV/c) but 
do increase rapidly from peripheral to central collisions. These scaling patterns have 
been interpreted as an indication that viscous damping for the higher harmonics follow 
the expected acoustic (or n 2 K) scaling [6]. Interestingly, these scaling patterns are also 
reflected in quark number (n q ) scaling of particle identified data. That is, v n /(n q ) n / 2 for 
different particle species, plotted as a function of transverse kinetic energy KEt, gives 
essentially a single curve as illustrated for f3/(n g ) 3//2 in the left panel of Fig. |3l The 
centrality dependence of t>3/(t>2) 3//2 and t>4/(f2) 2 also provide a valuable constraint for 
eccentricity models [6]; they favor a Glauber initial eccentricity model. 
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Figure 3. Illustration of quark number scaling for V3 (left). Comparisons of v n {^> n } 
vs. N part with theory as indicated (right). 



The right panels of Fig. [3] compare the centrality dependence of ^{^j an d ^{^s} 
with several hydrodynamical model calculations [91110], which demonstrates an essential 
constraint that the data provide. Panel (a) shows good agreement between data and 
calculations which employ MC-Glauber and MC-KLN initial eccentricities (£2) paired 
with viscosity values of 47r^ = 1 and 2, respectively The resulting rj/s uncertainty 
(~ 100%) reflect the model dependence of e 2 - Differences between the calculations 
become more apparent for the higher selection shown in panel (b), but the implied 
uncertainty for rj/s remains large. Panels (c) and (d) demonstrates the utility of ^{^3} 
as a constraint for disentangling the effects of viscosity and initial conditions; they 
indicate excellent agreement with the results from viscous hydrodynamics which employ 
Glauber initial eccentricities and Att'^ = 1, and rather poor agreement with calculations 
which employ MC-KLN initial conditions and Att^ = 2. The constraining power of 
f 3 {^/ 3 } stems from the fact that viscous corrections to v n scale as n 2 and the two 
eccentricity models give similar values for £3 [7]. Thus, the larger viscosity needed for 
agreement with the data with MC-KLN eccentricities in panels (a) and (b), leads to a 
significant under prediction of v 3 {^3} in panels (c) and (d). In summary, PHENIX v n 
measurements provide important constraints for robust extraction of rj/s. 
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